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I. INTRODUCTION 



Chiral effective field theory allows for a consistent calculation of forces between two, three and more nucleons (for 
recent reviews, see Refs. [lJ-Q ) . In this paper, we fill in the last missing part of the nuclear forces calculated at 
next-to-next-to-next-to-leading order (N 3 LO) in the chiral expansion, namely the short-range contributions and the 
leading relativistic corrections to the three-nucleon force (3NF). At this order, the 3NF is given by five topologies as 
shown in Ref. [4]. Three of these do not involve multi-nucleon operators and thus contribute at long range r ~ 1/M X 
and r ~ 1/(271/,,-), with M„ the pion mass. The corresponding momentum and coordinate-space representations are 
given in that paper (see also Q). Here, we work out the remaining terms corresponding to the two topologies involving 
four-nucleon operators (the short-range terms) as well as the relativistic (l/m) corrections to the leading (N 2 LO) one- 
and two-pion exchange contributions. Throughout, m denotes the nucleon mass. Special care has to be taken to 
calculate these corrections consistently with the ones contributing to the two-nucleon force. With the formalism given 
here, one is now in the position to perform calculations in three (or more) nucleon systems consistently at N 3 LO. It is 
important to stress that the 3NF at this order does not involve any new and unknown low-energy constants (LECs) - 
the full 3NF to this order thus depends just on two LECs, commonly called D and E, that parameterize the leading 
one-pion-contact term topology and the six-nucleon contact term at N 2 LO. Note that an exploratory study of the 
effects of the long-range parts of the 3NF at N 3 LO in the triton was recently presented in Q. 

Our manuscript is organized as follows: In Sec. [TT] we discuss the contributions of the one-pion-exchange-contact 
topology. Sec. [Till describes the calculation of the two-pion-exchange-contact topology while Sec. [IV] is devoted to 
the analysis of the leading relativistic corrections that also appear at N 3 LO. We end with a brief summary and 
outlook. For the sake of completeness, in appendix |A"1 we give the formal algebraic structure of the parts of the nuclear 
Hamiltonian which are relevant for our work. Finally, appendix [B] contains the coordinate-space representation of the 
obtained results. 
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FIG. 1: Diagrams contributing to the one-pion-exchange-contact 3NF. Solid and dashed lines represent nucleons and pions, 
respectively. Solid dots refer to the lowest-order (dimension one) vertices from the chiral Lagrangian. Diagrams which result 
from the interchange of the nucleon lines and/or application of the time reversal operation are not shown. 

II. THE ONE-PION-EXCHANGE-CONTACT TOPOLOGY 



We begin with the discussion of the one-pion-exchange contributions involving short-range contact interactions between 
two nucleons which are shown in Fig. [TJ It is important to keep in mind that these diagrams do not correspond to 
Feynman graphs. Rather, they represent schematically various contributions to the connected, irreducible part of the 
three-nucleon amplitude which gives rise to the three-nucleon force. Each diagram is to be understood as a set of all 
possible time-ordered-like graphs of the same topology (i.e. the same sequence of vertices). The precise meaning of 
the diagrams and their contributions to the nuclear Hamiltonian depend, strictly speaking, on the approach used to 
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extract the irreducible part of the amplitude. 1 In this work, we adopt the method of unitary transformation which 
was already successfully applied to the derivation of nuclear forces Ref. 04111 an d more recently electromagnetic 
nuclear exchange current and charge operators [H, [HI . It is also the same method that we used in Ref. Q to compute 
the long-range parts of the 3NF at N 3 LO. The resulting two- and many-body forces as well as the exchange currents 
are thus derived in the same theoretical approach from the effective chiral Lagrangian and are, per construction, 
consistent with each other. 

Before discussing the results for individual diagrams we first outline the way the actual calculation is carried out. 
Starting from the effective chiral Lagrangian for pions and nucleons, we first compute the corresponding Hamilton 
density using the canonical formalism, see [14| for more details. In order to arrive at nuclear forces (and current 
operators) valid in the low-energy region, i.e. well below the pion production threshold, the pion fields need to be 
integrated out. This is achieved through decoupling of the purely nucleonic subspace of the Fock space from the 
rest via a suitably chosen unitary transformation (UT) This step is carried out perturbatively utilizing the 

chiral expansion and making use of the standard chiral power counting based on naive dimensional analysis. We 
also exploit the freedom to choose the basis states in the nucleonic subspace of the Fock space and include a large 
number of additional unitary transformations which can be constructed at the given chiral order It turns out 
that the resulting unitary ambiguity of the nuclear Hamiltonian and exchange currents is strongly constrained by 
requiring that the resulting nuclear potentials and currents are renormalizable, see [Tl| for a detailed discussion. The 
explicit form of the (purely strong) unitary transformation compatible with the renormalizability requirement and the 
operators contributing to effective nuclear Hamiltonian which are needed in the present calculation can be found in 
appendix |A] see also Ref. . Once the unitary operator is determined at the desired order in the chiral expansion, 
the effective, purely nucleonic Hamiltonian can be computed straightforwardly. One ends up with terms given by a 
sequence of vertices and energy denominators which are similar to the ones emerging in time-ordered perturbation 
theory. However, the energy denominators and the coefficients in front of each terms generally differ from the ones 
arising in the context of time-ordered perturbation theory. We refer to [15| for more details on the method of unitary 
transformation. 

We now discuss the individual contributions to the three-nucleon force. Unless not stated otherwise, the expressions for 
the 3N potentials are to be understood as matrix elements with respect to the nucleon momenta while as operators with 
respect to spin and isospin quantum numbers. We begin with diagrams (l)-(6) in Fig. Q] with one self-energy /vertex- 
correction insertion on the second nucleon. The contribution of these diagrams to the 3NF depends on the unitary 
transformations given in Eqs. (3.24), (3.25) and (3.48) of Ref. Fixing the corresponding rotation angles 0:1,2,6 in 
the way compatible with renormalizability as described in that work leads to the following result for the corresponding 
3NF contribution 



g%C T 0W1 f dH 1 r 



DR 9a c t , , (<?i • qi)(a 3 ■ qi) 



lferJ? AJ ' Tl - T2 q\ + Ml ' (2 - 1} 

where u\ (tj) refer to the Pauli spin (isospin) matrices while q^ = pi ' — pi denote the momentum transfer of the 
nucleon i. Further, gA and F v stay for the nucleon axial-vector and pion decay coupling constant, respectively, while 
Ct is the low-energy constant accompanying the lowest-order (spin-dependent) two-nucleon contact interaction. We 
emphasize that the Wigner-symmetry-invariant two-nucleon contact interaction proportional to Cs yields a vanishing 
contribution which can be traced back to the fact that it commutes with the pion-nucleon vertex. Here and in what 
follows, we use dimensional regularization (DR) in the derivation of the 3NF. Further, we adopt the same notation as 
in Ref. Q and give results for a particular choice of nucleon labels (unless not stated otherwise) . The full expression 
for the 3NF results by taking into account all possible permutations of the nucleons 2 , i.e.: 

^3N U = ^3N + all permutations . (2.2) 

Next, the contribution of diagrams (7)-(9) in Fig. [T] with one self-energy/ vertex-correction insertion on the third 
nucleon is found to vanish completely. This feature does not depend on the structure of vertices appearing in these 
diagrams and emerges from the renormalizability constraints, see Ref. On the other hand, diagrams (10)-(15) 



1 We remind the reader that nuclear potentials do not correspond to observable quantities and are, in general, scheme dependent. 

2 For three nucleons there are altogether 6 permutations. 



FIG. 2: Examples of diagrams oc g\C% T which lead to finite shifts of the purely short-range 3NF at N 2 LO as explained in the 
text. For notation see Fig. [1] 



with pions being exchanged between the nucleons 1 and 2 and 2 and 3 do yield a nonvanishing contribution of the 
form 



V 



/d 3 l 1 r r -> -» i -> -o 



4F4 9 2 +M 2 
DR ffl^T (?1 ' <fl 
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M« 2 2 



2t1 • T 3 ((fl • a 3 ) - 7l • [r 2 X T 3 ] <fl • [(7 2 x <7 3 ] 



(2.3) 



We find that diagrams (16)-(18) involving a single self-energy/vertex-correction insertion on the nucleon 1 do not 
generate 3NF contributions. For graph (17) this feature is quite general and does not depend on the vertex structure. 
Irrcgardlcss of the choice for the additional unitary transformations, this diagram appears to be purely reducible at 
the order considered. Its contribution to the scattering amplitude is, therefore, properly accounted for by iterating 
the dynamical equation. The two self-energy corrections do yield non-vanishing irreducible contributions of opposite 
sign when considered separately. Since the leading nucleon self-energy correction is scalar, isoscalar and independent 
of the nucleon momentum, the contributions of these two diagrams sum up to zero. 

Similarly, diagrams (19)-(25) proportional to g\Cs,T do not produce any non-vanishing 3NF contributions. This is 
because the integrands entering the corresponding loop integrals are always odd functions of the loop momentum. This 
feature depends crucially on the renormalizability constraints discussed above which ensure that the pion-exchange 
between the nucleons 1 and 2 factorizes out in all irreducible contributions emerging from these diagrams. Thus, 
the pion loop integrals are always proportional to w ; -2 (which can be understood from the dimensional analysis) and 
involve a single power of the loop momentum I in the numerators emerging from the leading (derivative) pion-nucleon 
vertex ~ gA- 

Next, diagrams (26), (27) in Fig. [T] proportional to g\C\ T also lead to vanishing 3NF contributions for the choices of 
the additional UTs compatible with the renormalizability constraints. 3 It should be emphasized that there are also 
purely short-range contributions to the 3NF oc g\C\ T emerging from diagrams shown in Fig. [5] with NN contact 
interactions acting between different pairs of nucleons. While these graphs (presumably) generate non- vanishing 3NFs, 
their contributions are purely short-range and only provide a finite shift oc M v to the LEC E which accompanies the 
contact 3NF at N 2 LO. Since we anyway work with bare LEC E which needs to be refitted to experimental data at 
each order in the chiral expansion, there is no need to explicitly evaluate these contributions unless one is interested 
in the quark mass dependence of the 3NF and few-nucleon observables. 

Finally, the last two diagrams in Fig. []] clearly do not yield any contributions in a complete analogy with NLO one- 
pion-exchange-contact diagrams. Finally, there are also no contributions at N 3 LO from tree diagrams involving one 
insertion of the higher-order (^-vertices in the effective Lagrangian, see graphs (6) and (7) in Fig. [3] except for the 
relativistic corrections which will be considered in section IIVI As explained in Ref. [|[ , the contributions from these 
diagrams are suppressed by at least one power of Q/m where Q denotes a genuine soft scale. 



We are thus left with Eqs. (|2.1[) . (|2.3[) as the only nonvanishing contributions to the one-pion-exchange-contact 3NF 
topology. We now show that these terms cancel each other exactly if one takes into account the antisymmetric nature 



To avoid possible confusion, we emphasize once again that these diagrams correspond to all possible non-iterative time-ordered graphs 
wich nucleons treated as static sources. Non-static corrections are computed perturbatively as discussed in section ITVl 
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FIG. 3: Tree diagrams contributing to the two-pion-exchange and one-pion-exchange-contact topology of the 3NF at N 3 LO. 
The solid boxes denote insertions of either sub-subleading di-vertices from the effective pion-micleon Lagrangian or the leading 
1/m-corrections. For notation see Fig. [1] 

of few-nucleon states. In particular, we use the identities 

(V3CT3 + T 2 <T 2 2 ) = \ (T3CT3 + T 2 <T 2 - T 3 a 2 - Tia\ +T 2 X I-3 [<72 X CT3 ]*) 

\ / A23 4 

= B\ 

(t 2 X T 3 \S 2 X (73 f) = 2B l , 

V / A23 

(rs4 + r 2 4) A23 = -B', (2.4) 

where the superscript i refers to the cartesian component of the Pauli spin matrices and A23 denotes antisymmetriza- 
tion with respect to nucleons 2 and 3, which, for a momentum-independent operator X, can be written in the form 

(X) A23 .i(x-l±f^i±^x). (2.5) 

It is easy to sec that adding the contribution from interchanging the nucleons 2 and 3 to Eqs. (|2.1I) and (|2.3p and 
performing antisymmetrization with respect to these nucleons leads to a vanishing result. We, therefore, conclude 
that there are no one-pion-exchange-contact terms in the 3NF at N 3 LO. 



III. THE TWO-PION-EXCHANGE-CONTACT TOPOLOGY 



We now turn to the two-pion-exchange-contact diagrams shown in Fig. 3J Evaluating the matrix elements of the 
operators listed in Eq. (jA.ip for diagrams (l)-(7) in this figure we find the g^Cs.T-contribution to the two-pion- 
exchange-contact topology of the form 



V 2 „- 



g\C T f d 3 l 



m J (2tt 



Tl • T 2 



{ 4 2 + 24 ) (lliqi ■ °2 • #3 ) + 2qll 2 (a 2 ■ (7 3 ) q\{l ■ o 2 )(l ■ a 3 ) 



- qt(°2 ■ *s ) - l 2 (Qi ■ °2 ){qi ■ °z ) + l 2 (l • o 2 )(l ■ a 3 ) - l 4 {a 2 ■ a 3 )J + ^-j 
x ( - qf(qi ■ 02 <? 3 ) + q\{l- a 2 )(q\ ■ a 3 ) + l 2 (qi ■ a 2 )(f- a 3 ) - l 2 {l- a 2 ){q x ■ a 3 )) 
+ 2x 2 • r 3 ( - ir ^- + -^-f ) (q 2 l 2 (ai -a 2 ) + q 2 {l- o x )(f ■ o 2 ) + (qi ■ T)(qi ■ a x )(f ■ a 2 ) 

+ (9*1 cti ){q\ &2 ) (gl • l) 2 {&i ■ <? 2 ) - l 2 {qi ■ °i )(qi ■ <? 2 )) + 6( 4 1 2 + 2 X 4 

X ( - q\l 2 {Si ■ a 2 ) + ?i(f- CTr )(f- ff 2 ) - (<fi • f) (cfi • oi )(f ct 2 ) - (<fi • f)(f- ct x )(<fi • ct 2 ) 
+ (gl • 2 (^i • ^2) + J 2 (gi • <?i )(<fi • CT 2 ) 
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dr. 9aCt 



487rFi 



2xi • T2 (CT2 • <7 3 ) 



4M2 + 9 2 



+ 2{2Ml + ql)A{ qi ) 



(<fi ■ ai)(qi ■ a 2 ) - q x (oi ■ S 2 ) A(qi) } , 



where the loop function A(q) is defined according to 

A(q) = — arctan 
2q 



2AU 



(3.6) 



(3.7) 



Notice that we have exploited the fermionic nature of the nucleons in order to simplify the above expression and made 
use of the following identities: 



(r 2 + r 3 )a l 2 a J 3 



(r 2 ■ t 3 (ct 2 + a 3 



= --r(T 2 + T 3 )5ij(l-a 2 -a 3 ) 

A23 4 



A23 



= -3 (a 2 + <X 3 



A23 

3 3 
= (<?2 + ct 3 ) + -t 2 ■ t 3 (a 2 + o 3 ) , 



(3.8) 



where the ellipses in the first line refer to terms which are antisymmetric with respect to the interchange of the indices 
i, j and therefore do not contribute to the final result. 

We now turn to diagrams (8)-(12) in Fig. |4] whose contributions are proportional to g\Cs,T- We find that graph (12) 
does not contain any irreducible pieces, while the contributions of diagrams (8)-(ll) reads: 



V 2 „. 



9 2 A 



d 3 l 



32F* J (2tt) £ 

+ 4C T Tl ■ T 2 



,2 , ,2 



i(C s + C T ) [ti x t 2 \ ■ t 3 - (qx ■ a 2 )(l ■ a 3 ) + (I ■ a 2 )(q\ ■ a 3 ) 



q\(o 2 -a 3 ) + (qi- <j 2 ){q\ ■ a 3 ) + 1 2 {S 2 ■ a 3 ) - (I ■ a 2 ){l ■ a 3 ) - 2i[qx x I } ■ a 3 .(3.9) 



It is easy to see that all terms involving the imaginary unit number i are vanishing. Using dimensional regularization 
and performing antisymmetrization with respect to the nucleons 2 and 3 we arrive at the following final result for the 
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g\Cs,T contribution to V^-c 



DR 9 2 A C T 



247rFi 



Ti • t 2 (a 2 ■ a 3 ) M n + {2Ml + ql)A{ qi ) 



(3.10) 



Finally, it is easy to see that the last diagram in Fig. @] does not contribute to 3NF as the corresponding Feynman 
graph involves at this order in the chiral expansion only reducible pieces. To summarize, the complete contribution 
of the two-pion-exchange-contact topology of the 3NF at N 3 LO is given by Eqs. Q3.6P and (|3.10p using the convention 
of Eq. (IQl . 



IV. LEADING RELATIVISTIC CORRECTIONS 



The leading relativistic corrections - that is the corrections of the operators in 1/m - to the 3NF provide further 
contributions to the two-pion-exchange and one-pion-exchange-contact topologies. They emerge from two different 
sources. We remind the reader that the formally leading 3NF generated by the tree-level two-pion-exchange 4 and 
one-pion-exchange contact diagrams vanishes at NLO. Stated differently, the resulting contributions to the scattering 
amplitude are either purely reducible (at that order) or shifted to N 3 LO due to the suppression by one power of Q/m 
caused by the time derivative entering the Weinberg- Tomozawa vertex [H|,[l3]- The first kind of relativistic corrections 
emerges from taking into account retardation effects in NLO diagrams, see graphs (1), (2) and (6) in Fig. [3] As will 
be shown below, contrary to the case when the nucleons are treated as static sources, these diagrams do produce 
non- vanishing 3NFs. Secondly, one also needs to take into account the 1/m-corrections to the leading ttNN and 
ttttNN vertices which leads to graphs (3), (4), (5) and (7) in that figure. Notice that there are no 1/m-corrections to 
the leading two-nucleon contact interactions. 

We begin with the retardation corrections to the two-pion exchange 3NF corresponding to diagram (1) in Fi g. EQ 
When evaluating its contribution to the potential, one needs, in addition to the terms listed in Eq. (A.l) of Ref. [llj j . 
to take into account effects induced by the additional unitary transformation driven by the operator S$ in Eq. (4.24) 
of Ref. [13J which is parametrized in that work by a constant /3g ■ The explicit form of the operators to be evaluated 
is given in Eq. (|A.3[) . The resulting 3NF contribution has the form 



Vn 



2tt, 1/r, 



(<7l ■9l)(o ! 3 -<f 3 ) 



(1 - 2/3 8 )(-ri • t 3 (qi ■ q 3 ) 2 + [ti x r 2 ] • 1-3 [qi x q 3 ] ■ a 2 (qi • $3) 



32mF4(^ + M2)2( g | + M2) 
2i(i-i • r 3 [ql x q 3 ] ■ a 2 - [n x r 2 ] • t 3 (q\ ■ <f 3 ) ) ( (1 - 2j3 8 ) (<fi • k 2 ) + (1 + 2/3 8 ) (<fi • k x ) 



(4.11) 



where ki = (pi +pi')/2. Notice that the unitary transformation considered above also affects the form of the 1/m 2 - 
corrections to the one-pion exchange and 1/m-corrections to the two-pion exchange two-nucleon potentials at N 3 LO. 
To be consistent with the N 3 LO potential of Ref. 10;], one has to choose (3$ = 1/4. The contribution of diagram (2) 
in Fig. [3] is not affected by the above UT and has the form 



V: 



2tt, 1/r 



9\ 



(gi ■ qi Kgg • 93 ) 



[t x x t 2 ] ■ r 3 (q 3 ■ k 3 - qi ■ ki) . 



(4.12) 



Finally, retardation corrections to the one-pion-exchange-contact topology from diagram (6) in Fig. [3] have the form 



cont, 1/m 



9l 



<7i • qi 



8mFZ(ql + M*y 



TfT 



2 {(1 - 2/3 8 ) (qi ■ q 3 ) [C s {qi ■ a 2 ) + C T {q x ■ a 3 )] + 2iC 



(1 - 2^ 8 ) (q x ■ k 2 ) + (1 + 2^ 8 ) (gi • fc x )] } 



rqi ■ [02 x <j 3 \ 

(4.13) 



4 There are two such diagrams: the one proportional to g^ A and the one involving a Weinberg- Tomozawa vertex which is proportional to 
9 2 A- 
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We now turn to contributions emerging from relativistic corrections to vertices in the pion-nucleon Hamilton den- 
sity. Consider first the contributions of diagrams (3), (4) and (7) in Fig. [3J The relevant terms in the effective 
Lagrangian/Hamiltonian have the form 5 

C wN = -H wN = -M-rfrff ■ VttJV - i^ r NWa • (V - V) N ■ -k , (4.14) 

where N and tt refer to the nucleon and pion fields, respectively. Similarly to the retardation corrections considered 
above, it is possible to construct an additional UT with the generator proportional to g\/m which affects the contri- 
butions of diagrams (3) and (7) when acting on the one-pion exchange and the leading contact potentials, respectively. 
We adopt here the notation of Ref. [13J , where this unitary operator is parameterized in terms of another arbitrary 
parameter /3 9 , see Eq. (4.24) of that work. Again, the UT also affects the 1/m 2 - and 1/m-corrections to the one-pion 
and two-pion exchange potential at N 3 LO, respectively, so that the parameter (3g needs to be chosen consistently. 
Using the explicit expressions for the operator structure of the unitarily transformed Hamilton operator in Eq. (A. 4) 
of Ref. [HI and adding terms induced by the above mentioned UT we obtain the following 3NF contribution from 
diagram (3): 

x [qi x q 3 ] ■ a 2 ] + 2i[n x r 2 ] • r 3 [ - (2/3 9 - 1) (a 3 • <? 3 ) («i • h) + (2/3 9 + 1) (<r 3 • fc 3 )(<fi • q 3 )] \ . (4.15) 



Next, it is easy to see that diagram (4) does not generate any 3NF at the considered order in the chiral expansion. 
Similarly to the corresponding NLO diagram (same topology constructed from the lowest-order vertices), there is 
suppression by, in this case, two factors of Q jm due to the time derivatives entering the Weinberg- Tomozawa vertex 
and the 1 / m-correction to the ttNN vertex in Eq. (|4.14l) . Consequently, within the power counting scheme adopted 
in the present work, this diagram contributes to the 3NF only at order N 5 LO. Finally, the contribution of diagram 
(7) in Fig. [3] is clearly also affected by the abovementioned unitary transformation and is given by 

q\C t 1 f 

ViTr-cont, i/m = ~~ Q — -E^^TT T7I T i ' x 2 { 2i(2/3 9 + 1) (fci • a ± ) q\ ■ [a 2 x <r 3 ] - (2/3 9 - 1) (<fi • ffi) (q 3 ■ a 3 ) (4.16) 
8m_ri at + Mi L 



8mF% q\ + Ml 

2i(2j3 9 - 1) (q\ • 3 X ) k 2 ■ [a 2 x a 3 ]\ + 9 } G * 2 1 T t • r 2 (2^ 9 - 1) {q\ ■ a x ) (q 3 ■ a 2 ) . 

J 8ml 1 * q{ + M£ 



Finally, consider the contribution of diagram (5). Since this graph does not induce any reducible pieces, the cor- 
responding 3NF can be identified with the amplitude and computed using the Feynman-diagram technique. The 
Feynman rule for the 1/m-correction to the ttttNN vertex can be found e.g. in Ref. [13]. In our notation the relevant 
terms 6 have the form 



, mF f~T- ip + f) • (.I + q 2 ) - ^rf^ X q 2 ] ■ ff, (4.17) 

where q±, a (q 2 , b) denote the momentum and isospin quantum number of the incoming (outgoing) pion. We then 
obtain for the contribution of diagram (5): 

VW, 1/m = I 2 S l\flfl'ft 2 , [Tl X T 3 ] ■ T 2 ([q, X q 3 ] ■ a 2 + l k 2 ■ {q 3 - q,)) . (4.18) 

32mF| (qi +M£){qi + M%) V / 

To summarize, our results for the relativistic (1/m) corrections to the 3NF at N 3 LO are given by Eqs. (|4.11[) . (|4.12p . 
(I4.13p . (|4.15p . (|4.16p and (|4. 18|) . Adding up these expressions, the two-pion exchange contribution can be written in 
the form 



Y _ 9a 1 

V2-K, 1/r: 



32mF* {ql+Ml){ql+MD 



\T\ ■ t 3 Ai 23 + [n x t 2 ] ■ t 3 Si 23 ^ + 5 permutations , (4-19) 



5 Different conventions for the sign of the axial vector coupling constant gA are used in the literature. For example, gA in Ref. [l8| 
corresponds to —gA m the present work. This difference does not affect the expressions for nuclear forces and currents where the LEC 
gA enters quadratically. It is important, however, to use the same convention for both terms in Eq. I|4.14ll . 

6 We do not show terms involving zeroth components of the nucleon momenta since they do not contribute to the 3NF at N 3 LO. 
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where 

A 12 3 = (#1 ■ qi)(&3 ■ 93) «123 + (<?1 ' <Il)(&3 ■ fa) C123 , 

-B123 = (Si • q\)(S 3 ■ 93)^123 + (<?i • ?i)(o*3 • fa)di23 , 
and the functions ai 23l 6123, C123 and c£i23 are given by: 

(1 - 2fo)(qi ■ qs) 2 - 2i[q\ x q 3 ] ■ a 2 ((1 - 2fa)(q l ■ k 2 ) + (1 + 2£ 8 )(<f 1 • 

2z[gi x k 2 ] ■ S 2 
l)[<fi x 33] • <?2 , 

(1 - 2/3 g j |>/, X ,/;, 



(4.20) 



Q-123 



ql 



Ml 



Cl23 
&123 



1) 



3? 



3? 



M 2 



3) Wl X 33] • <?2 (qi ■ q 3 ) + 2i{q x ■ q 3 ) ((1 - 2 ) S 8 )(<fl ■ k 2 ) + (1 + 2/3 8 )(gi • &)) 



^123 — 



+ 2ig 3 • (fc 3 - fe 2 ) + 2ig 2 A (2p 9 - • fc 2 ) - [ft x g 3 ] 
-2^(209 + !)(& -q 3 ). 



<?2 , 



The relativistic corrections to the one-pion-exchange-contact topology have the form 

1 



cont, 1/m 



8mF 2 q\ + Ml 



Tl • T 2 ( (CTi • (ft) /i 23 + (ffi • fci) gi 



23 



5 permutations . 



(4.21) 



(4.22) 



where 
/123 

5123 



1 



ft 2 



M2 
x (gl ■ fci) 
-2iC T (20 9 



( I ~ 2. k )( 7, • ./V) ( r\-l -7, • ^ ) + f , ! <7i ) + 2tC T fi • [* 2 * <r 3 ] ((1 - 2&)(gi • fe 2 ) + (1 + 2ft) 

2iC T (2p 9 - 1% ■ [S 2 x ^ 3 ] , 



+ (2/3 9 - 1) ]C s (q 3 ■ S 2 ) + C T (q 3 ■ ffa) 
1)9*1 • [S 2 x S 3 ] . 



(4.23) 



The above expressions are written in a general form including the dependence on the constants /3s. 9 which parametrize 
the unitary ambiguity of these potentials. The 1 /m 2 -corrections (1/m-corrections) to the one-pion exchange (two- 
pion-exchange) two-nucleon potential at N 3 LO are also affected by this unitary ambiguity. A comprehensive discussion 
on this issue and on the relation between different forms of relativistic extensions of the Schrodinger equation and non- 
static terms in the two-nucleon potential can be found in Ref. [20]. To be consistent with the two-nucleon potential 
of Ref. Oil , one needs to choose 



09=0. 



(4.24) 



We emphasize that in order to completely take into account the relativistic corrections at N 3 LO in few-body calcu- 
lations, one needs to use the Schrodinger equation with the relativistic expression for the nucleon kinetic energy, the 
proper forms of the 1 / m 2 -corrections to the l7r-exchange NN potential as well as 1/ m-corrections to the NN and 3N 
forces. In addition, one has to take into account boost corrections to the leading-order two-nucleon potential (i.e. the 
P- dependence in the i7r-exchange potential and the leading contact interactions, see Ref. (2l|). 

While 1/m-corrections to the one-pion-exchange-contact 3NF given in Eq. (|4.22[) have been never considered before, 
the leading relativistic corrections to the two-pion exchange 3NF were studied by Friar and Coon in Ref. [22j, who, 
however, only provide expressions in coordinate space. In order to facilitate a comparison between our results and the 
ones by Friar and Coon, we give in appendix iBlthe coordinate-space representation of the obtained 3NF contributions. 
Our profile functions are related to the one used in Ref. [22| according to 



47T 

Ui(M„r) = j^h (r) 



V P U 2 (M^r) = -AirM„rh (r) 



(4.25) 



where the function h is defined in Eq. (24c) of that work. Further, f 2 (M) in Ref. [22] correspond to (g^/(2F 7r )) 2 
(m) in our notation. Friar and Coon also discuss the unitary ambiguity of the resulting 3N potentials associated with 
the UTs oc 1/m. They parametrized it in terms of the arbitrary parameters /i and v, see also Ref. 20]. The relations 
between their /x, v and our /3s, 9 are given by 



\i = Af3 9 + 1 , 



2/3 8 



(4.26) 
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Taking into account the typographical error in [22j as pointed out in [5] of Ref. [20(, we observe that our results 
for the 1/m-corrections to the two-pion exchange 3NF in Eqs. (|B.3p - (|B.5p agree with the ones in Eqs. (23a), (23b), 
(23c), (33a), (33b), (33c) and (33d) of Ref. [22[ except for the retardation corrections in a\2z and 6123 emerging from 
diagram (1) in Fig. [3] More precisely, we reproduce the results of Friar and Coon if we make the replacement 

((1 - 2 J 8 8 )(Si • k 2 ) + (1 + 2A,)(& • fc x )) — ► ((1 - 2 J 8 8 )(fi • k 2 ) + (1 + 2 / S 8 )(£ • hj) + S (4.27) 

in 0123 and 6123 in Eq. (14.21)) with 8 given by 

S = (1 - 2&)gl • (k 2 - h) + (1 - 2fo)q 3 ■ (k 2 - k 3 ) = i (1 - 2/3 8 ) (p? + pj + p\ - p[ 2 - p' 2 2 - p' 2 ) , (4.28) 

where we made use of the identity qi + q 2 + q 3 =0. This implies that the difference between our results appears 
only when the nucleons are off-the-energy shell. We conjecture that this difference originates from the additional UTs 
exp(aiSi + a 2 S 2 ) with ai.2 denoting the transformation angles and S\,2 the generators of the two-pion-exchange range 
whose explicit form is given in Eq. (3.25) of Ref. [111 ]. The corresponding unitary ambiguity is not explored by Friar 
and Coon. It is easy to see that these UTs induce, among others, terms in the 3NF oc [E'kin, (ceiSi + 0:2 S2)] which 
have the same structure as the <5-terms discussed above. As explained in Ref. [TTj] , the angles ai :2 have to be chosen 
in a specific way, see Eq. (3.31) of that work, in order to maintain renormalizability of the leading loop corrections to 
the 3NF. Our findings therefore indicate that the particular choice made in Ref. [22| for ot\. 2 is not compatible with 
the abovementioned renormalizability constraint. 



V. SUMMARY AND CONCLUSIONS 



In this work, we have derived the short-range part and the relativistic (1/m) corrections to the 3NF in chiral effective 
field theory at N 3 LO. Combined with the long-range parts already given in Ref. this completes the calculation of 
two-, three- and four-nucleon forces at this order in the formulation without explicit A(1232) degrees of freedom. The 
short-range parts considered here consists formally of two different topologies. Remarkably, there is no contribution 
from the one-pion-exchange-contact topology at N 3 LO as we have demonstrated in Sec. [TTJ The total contribution 
from the two-pion-exchange-contact topology is given in Eqs. p.6[) , (|3.10j) supplemented with Eq. (|2.2p . Further, there 
are relativistic {1/m) corrections to the leading 3NF at N 2 LO, which are worked out in Sec. IIVI The corresponding 
terms of the 3NF are summarized in Eqs. (14. 19|) and (I4.22p . We also provide the coordinate-space representation 
of the obtained results in appendix [B] and compare our findings for the long-range 1 / m-corrections with the earlier 
calculation by Friar and Coon [22j . 

We stress again that at N 3 LO, the 3NF is free of unknown LECs, so the complete 3NF consisting of terms at N 2 LO 
and N 3 LO contains altogether only two LECs D and E which need to be determined from few-nucleon data. With 
the results presented here, one is now in the position to analyze in detail the many data in few-nucleon and heavier 
systems based on consistent and precise two- and three-nucleon forces. In particular, it will be interesting to find 
out whether the remaining discrepancies in the three- and four-nucleon systems will be overcome employing the force 
derived in Ref. [J] and here. Work along these lines is in progress. 
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Appendix A: Formal algebraic structure of the 3NF corrections 



In this appendix we list the formal operator structure of the various N 3 LO contributions to the nuclear Hamiltonian 
relevant for the present calculations. A detailed discussion on these terms can be found in Re f. fill) - We include 
contributions induced by the additional UTs considered in that work and in Eq. (4.23) of Ref. [1J|. Except for etc, 
and /3s, 9> the corresponding transformation angles are fixed by the renormalizability constraints as discussed in (Ilj |. 
In all applications to nuclear forces and currents considered so far, the dependence on 015 drops in the final results. 



• terms oc g\Cs,T'- 



V = r\ 



1 f/W Al ffW r, f/W Al Al f/W 1 TT& Al f/W r, f/W Al Al 
- 2 H 2i ~e~ H 2i 7 l tl 2i ~e~ H 40 ~E2 H 2i ~ 2 21 ~E~ 21 1 21 E 2 40 E~ ' 

2 21 El 21 ' 21 E„ 40 E„ - 
A 1 



21 



J 21 



1 ttW r, h-W A A h-C 1 ) A * h-M 



*21 



2 40 / 21 ^ 21 p^ 21 E 2 21 
1 H"( 2 ) r. WW Al rr(l) ^ rr(l) ^ rr(l) 

-H i0 r)H 21 ^2 #21 e~ m 21 e~ H ^ 



(1) Al rr(l) 

91 J^91 



1 H^KfJ^ Al H^KfJ^ Al 
2^40 ^#21 £"-"21 ?7#21 



A 2 



1 ,, I A rr(l) f/ - 
9 21 "p~#21 ~p~#40 TT 

^ fi^r 



3 r, TT^ Al f/W « fT^ Al 

-H i0 T)H 21 -g2-H 2 l ^#21 ^2^21 

^" 4 n * #21 



-iff 



7T 7T 

rr(l) A 1 „(1) A 2 rr(X) A 1 „(2) A 1 , 

H 21 TT H 21 TT H 21 TT H 40 ~p~ < 

f-J-K 

(2) A 2 (i) A 1 (i) (1) A 1 (i) (i) A 1 (i) ( 2 

40 £~ 21 ~E~ 21 + 21 ~E~ 21 1 21 £?3 21 ^ 40 

4 21 £ T 21 £ T 21 £ 2 21 ' 40 



-(2) 
40 



1 Al h-W A2 wW Al r-W w( 2 
4^21 ^21 ^21 ^#21 VH 40 



-h.c. 



(A.l) 



terms oc g\Cs,T' 



V = n 



rr(2) 
"#22 



+ 



o n tt&KhW Al X * 
- 2a 5 H i0 r]U 21 -£T H 2i ~gi 

' X j. n \ » ff< 2 ) A2 ffW Al 

- +Q-5 1 Yi 40 ?? /i 2 2 ~Jjj~ M 21 ^2^21 



- + Q'. 



, ^ tt^KhW Al h-< 2 ) Al W^ 1 ) 

: 5 I #40 77 #21 ^~ W 22 ~~jP2 H 21 



(i) A 1 
21 £? H 



(i) A 2 „(2) 
21 p 



3 22 



-(2) „(1) A 1 (2 ) A 1 (1) 

40 '/"21 g2 22 p "21 

\ 1 - ' 

21 Tp 



f1 MVn, } ff (2) r, ?T (2) A2 Al 

(l + 2a 5 J-H 40 77^22 -g2~ H 2i ~jf H 2i 

TP 4U n 21 tp 11 



7T — " 

/HI) A 1 rr(l) A 2 (2) A 2 ( 2 ) _ (1) A 1 ( 2 ) A 1 ( 2 ) A 1 (1) _ (!) A 1 ( 2 ) A 1 (!) A 2 
"21 p -"21 p -"40 p "22 "21 p -"22 p "40 p "21 "21 p -"40 p "21 p 



V 



(A.2) 
-h.c. 



• retardation corrections oc g\/m: 



V = r/ 



1 4- fl ^ H-f 1 ) Al H-f 1 ) r, PfW Al H-^ 1 ) r, Iff 1 ) A 1 rr(l) _ A 1 (2) . 1 „(1) 

2+P8 I #21 ^^21 ^#21 ^3^21 ? /#20 _ #21 ^^21 ^#21 ^3^20 A #21 

^ X ^ f/W Al H-t 1 ) « f/( 2 ) « f/^ 1 ) Al f/^ 1 ) 3 wC 1 ) Al H-^ 1 ) ^ TT-^ 1 ' Al H-^ 1 ) « 77-( 2 ) 

^2 ~ / 21 £~ 21 ^ 20 V 21 £3 21 _ 4 21 21 £~ 21 ^2 21 ^ 20 

rjTX ^ rjTX -CJ-K 

.(1) A 1 (!) A 2 ( 2) 2 rr(l) A 1 (1) 3 (1 ) A 1 (1) (1 ) A 1 (l) ( 

>: —■//.,, — r ././., n A # 21 -gT#21 + g #21 ^2^21 ^#21 ^2^21 ^#2 



r(2) 
J 20 
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J 21 



~~ 2 'e 2 21 ^ 21 'e 2 20 1 

_ I iT« *1 — »? #1? 

4 21 £ 2 21 £L 21 £ 2 21 ' 20 



8 ~EP- 20 'E 2 21 

77 + h.c. , 



(A.3) 



retardation corrections oc g 2 A /m: 



rr(l) Al rr(l) ^ rr(2) x 2 rr(2) 
-"21 £i "21 g2 20 A 1:1 22 



(l + 2a 5)jff «A.if«^. 

1 C1 9™ X Al Al rj(2) rr(l) Al rr(2) A * rr(2) ,1 „(1) 

-(1 - 2a 5 j W 21 ^#22 £2#21 7 /#20 ~ #21 £^#22 E 2 H 20 A #21 

l(l + 2a 5 )H£^H£^H$ V H$ + i(l + 2a 5 ) !#>|^>|-fl£> „ ^ 



„(1) A 1 ( 2 ) x (1) A 2 (2 ) (2) A 2 (1 ) A 1 (!) (2' 

_ #21 £2 '-"20 A #21 £f#22 - <*a 5 H 22 ^"#21 ^#21 V #20 



A 2 

9n , N o-(2) A2 „(1) A 1 

-(l-2a 5 )H 22 -^H 21 — 



#21 7 7#20 



-(2) 



Al H^Kh^Kh^ 
- j '""21 £3^21 ?7#20 ^#20 



+ -"21 £3^20 A #21 ?7#20 ^ 2 21 ^3 20 20 21 



-(2) 
20 

h.c. , 



(A.4) 



retardation corrections oc g\Cs,T/in 



V = T] 



1 Al WW Al tfW r. rr(l) Al rr(2) Al rr(2) x l rr(l) , 1 rr(l) Al rr(2) Al „(1) „(2) 

2 #21 ^#40 ^2 #21 *7#20 _ #21 ^~#40 ^2^20 A #21 + 7J ^21 ^#40 ^#21 ? 7 #20 

f 1 fl ^ WW Al H"W r, FfW *, f/( 2 ) 4- f/C 1 ) Al fr(2) \1 rr(l)„ o-(2) 

I 2 8 / 21 21 ' 20 ^ 40 21 EP 20 21 ^ 40 



f 1 4- R \ TjW „ rr(l) Al tt(1) „ rr(2) 

~~ V 2 / 40 ^ 21 20 



r/ + h.c. , 



(A.5) 



• terms involving 1 /m-corrections to the <ju-vertex, cx g\ jm 

V - J f 1 J-R \ H-W Al H^KfJ^ Al ff( 3 ' '' ' 

v - V I 2 J ~E~ 21 ?1 ~ 



(1) Al rr(l) A2 rr(l) Al rr(3) _ rr(l) Al rr(l) ^ „(3) A 1 
21 77 -"21 77 -"21 TTi -"21 -"21 77 -"21 77 -"21 77 -"21 

&Tt &1T J^TT 



Al n (1 Kn^ Al 

#21 ^#21 T^"#21 



2 #21 ^21 ^#21 ^#21 



2 21 E 2 21 / 21 E 21 



(A.6) 



• terms involving 1/m-corrections to the c^-vertex, cx g\C s ,t I tn'- 

' (1) Al 77(2) A 1 (3 ) fl -\ (i) A 1 ( 3 ) (2) / 1 

#21 ^#40 ^"^21 - I 2 9 / 21 E 2 21 40 ~ \2 



V = ?? 



"^9 



rr(3) Al 77(1) „(2) 
#21 "^2"#21 ^#40 



r/+h.c. . (A. 7) 



Here and in what follows, we adopt the notation of Refs. [3. [Ill - IT3j . In particular, the subscripts a and b in refer to 
the number of the nucleon and pion fields, respectively, while the superscript k gives the inverse mass dimension of the 
corresponding coupling constant , see Ref. \m for more details. The chiral order associated with a given contribution 



For 1 /m-corrections, K; corresponds to the inverse power of coupling constants plus twice the power of m 1 . In particular, k = 2 for 
the nucleon kinetic energy term H20- 
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can easily be read off by adding together the dimensions k of . More precisely, it is given by J2i K i ~ 2. 111 the 
above equations, rj (A) denote projection operators onto the purely nucleonic (the remaining) part of the Fock space 
satisfying r/ 2 = rj, A 2 = A, 77A = Xrj = and A + 77 = 1. The superscript i of A 1 refers to the number of pions in 
the corresponding intermediate state. Further, E n denotes the total energy of the pions in the corresponding state, 



-Ett = Xli v h 2 + with ij the corresponding pion momenta. 



Appendix B: Coordinate-space representation 

In this appendix we give the coordinate-s pac e representation of the obtained 3NF contributions. Here and in what 
follows, we adopt the notation of Refs. [3, 11] and use the dimension- less profile functions Ui 2 , W and Wi, 2 defined 
for a general form of a local regulator function F A according to 

J (2ny q z + M£ x 



M 2 J (2tt) 3 vyy 2x 2 ' 

W2{x) = I Wf e ' m q2A{q) FA(q) = ~ v ' Wl{x) ■ (R1) 

For the two-pion-exchange-contact terms in Eqs. (|3.6l3.10j) we find: 

VW-cout = 19^1 { 2t i • T2 ^ 2 • ^) H 3 ^ - l ) W M - 2g 2 A U 1 (2x 12 ) + (2g 2 A - l)(2W 1 (x 12 ) + W 2 (x 12 )) 



-9g 



(ai ■ Vi 2 )(ct2 • Vi 2 )VFi(xi 2 ) + ((?!• a 2 )W 2 (xi 2 ) \w(x 32 ) + 5 permutations . (B.2) 



Next, the coordinate-space representation of the relativistic corrections to the two-pion exchange 3NF emerges from 
taking the Fourier transform of Eq. (I4.19P : 

g 2 ( \ 

V 2 „, i/ m = 1 Q24^ 717r 2^4 ( Tl - T 3^i23 + [ti x t 2 ] ■ t 3 Bi 23 ) + 5 permutations , (B.3) 



where 



^123 = (<?X • Vl2)(?3 ■ V32) 5l23 + I (^3 ■ j^-J, • V12) 6123} , 

-B123 = (cti • Vi 2 )(o ; 3 • V 32 )6i 2 3 + • (^i ' ^12)^123} j (B.4) 

and the functions 0123, &123, C123 and di 23 are given by: 

ai23 = g A (l-2p s )(V 12 - V 32 ) 2 C/ 2 (a;i2)C/i(x32) 

+ g 2 A {(l - 2ft)-jj|- + (1 + 2 A0|^, Vxa [V12 x V 32 ] • S 2 U 2 (x 12 ) Ui(x 32 )} - 2g A (2p 9 - 1)V 2 2 C^ifxia) U 1 (x 32 ) 

+ 2 5 i(2^ 9 -l){^, [^xV^t/iMC/i^)}, 
5i23 = -2<7i(2/3 9 + l)[Vi2 x V 32 ] -a 2 Ux{xi 2 ) t/i(x 32 ), 
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6i23 = g 2 A (l-2^ 8 )[V 12 x V 32 ] -a 2 (Vi 2 - V 3 2)t/ 2 (2:i2)C/i(x32) 

- 5 2| (1 _ 2( g 8) A + (i + 2y 8 8 )A j V 12 (V 12 -V 32 ) £/ 2 (a; 12 ) U 1 (x 32 )}-2{^--^, V32 t/i(*i 2 ) C/i^)} 

- 2. 9 ^ (2 A, - V12 (7i(x 12 ) C/i(x 32 )} - 2[Vi2 x V 32 ] • a 2 Ut(x 12 ) Ux(x 32 ) , 

rf 123 = 2 ff i(2A> + l)(Vl2-V 3 2)tfl(zi2)tfl(Z32), (B.5) 

where iy = M^r^ , Xij = \xij\ and ry = — is the distance between the nucleons i and j. Further, the Vy act on 
(i.e. are dimension-less) and refer to the momentum operator of the nucleon i. Curly brackets denote, as usual, 
the anti-commutator of two operators, {X, Y} = XY + YX, and {X, Y} = X ■ Y + Y ■ X. 

Finally, the 1 /m-corrections to the one-pion-exchange-contact topology have the form 

^lTT-cont, i/m = ^~f^ Tl ' T 2 (Vi • V i2 ) /i 23 + |?i ■ jj-, ffi 23 j) + 5 permutations , (B.6) 
with the functions /i 23 and gi 23 given by: 

/ 123 = (l-2^ 8 )(V 12 -V 32 )(Cs(V 12 -a 2 ) + CT(V 12 - ( 7 3 ))c/ 2 ( a : 12 )^(x 32 ) 

- C T {(1 - 2/3s)^ + (1 + 2/3 8 )^, Vi 2 V12 • [<J 2 x a 3 ] U 2 (x 12 ) W(x 32 )} 

- 2(2^ 9 - l)(c s (V 32 • d 2 ) + C T (V 32 ■ a 3 )) (7x(x 12 ) W(x 32 ) + 2C T {2% - p2 x ?s] ^12)^32)} , 
5123 = -2C T (2/3 9 + l)Vi 2 • [S 2 x <j 3 ] U x {x 12 ) W{x 32 ) . (B.7) 
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